Dissipative dynamics of two-qubit system: four-level lasing 
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The dissipative dynamics of a two-qubit system is studied theoretically. We make use of the Bloch-Redfield 
formalism which explicitly includes the parameter-dependent relaxation rates. We consider the case of two flux 
qubits, when the controlling parameters are the partial magnetic fluxes through the qubits' loops. The strong 
dependence of the inter-level relaxation rates on the controlling magnetic fluxes is demonstrated for the realistic 
system. This allows us to propose several mechanisms for lasing in this four-level system. 



I. INTRODUCTION 

Recently considerable progress has been made in studying 
Josephson-junctions-based superconducting circuits, which 
can behave as effectively few-level quantum systems. 1 When 
the dynamics of the system can be described in terms of two 
levels only, this circuit is called a qubit. Demonstrations of 
the energy level quantization and the quantum coherence pro- 
vide the basis for both possible practical applications and for 
studying fundamental quantum phenomena in systems involv- 
ing qubits. Important distinctions of these multi-level artifi- 
cial quantum systems from their microscopic counterparts are 
high level of controllability and unavoidable coupling to the 
dissipative environment. 

Multi-level systems with solid-state qubits may be realized 
in different ways. First, the devices used for qubits in reality 
are themselves multi-level systems with the lowest two levels 
used to form a qubit. For some recent study of multi-level 
superconducting devices see Ref. 2, Then, a qubit can be cou- 
pled to another quantum system, e.g. a quantum resonator. 3 
Such a composite system is also described by a multi-level 
structure. As a particular case of coupling with other systems, 
the multi-qubit system is of particular interest (see e.g. Ref. 

Operations with the multi-level systems can be described 
with level-population dynamics. In particular, population in- 
version was proposed for cooling and lasing with supercon- 
ducting qubits. 5,6 However, most of the previous propositions 
were related to three-level systems, while for practical pur- 
poses four-level systems are often more advantageous. 7 

The natural candidate for the solid-state four-level system 
is the system of two coupled qubits. The purpose of this paper 
is the theoretical study of mechanisms of population inversion 
and lasing, as a result of the pumping and relaxation processes 
in the system. We will start in the next Section by demon- 
strating the controllable energy level structure of the system. 
Our calculations are done for the parameters of the realistic 
two-flux-qubit system studied in Ref. 8. To describe the dy- 
namics of the system we will present the Bloch-Redfield for- 
malism in Sec. III. The key feature of the system is the strong 
dependence of the relaxation rates on the controlling param- 
eters. Then solving the master equation in Sec. IV we will 
demonstrate several mechanisms for creating the population 
inversion in our four-level system. We will demonstrate fur- 
ther that applying additional driving induces transitions be- 
tween the operating states resulting in stimulated emission. 



We summarize our theoretical results in Sec. V. and, based on 
our calculations, we then discuss the experimental feasibility 
of the two-qubit lasing. 



H. MODEL HAMILTONIAN AND EIGENSTATES OF THE 
TWO-QUBIT SYSTEM 

The main object of our study is a system of two coupled 
qubits. And altough our analysis bears general character, for 
concreteness we consider superconducting flux qubits, see 
Fig. Q] A flux qubit, which is a superconducting ring with 
three Josephson junctions, can be controlled by constant ("JMc) 
and alternating ($ ac sinwi) external magnetic fluxes. Each of 
the two qubits can be considered as a two-level system with 
the Hamiltonian in the pseudospin notation^ 

Hll^—e^- 1 -^, (1) 

where Aj is the tunnelling amplitude, are the Pauli ma- 
trices in the basis {|t), |t)} of the current operator in the z-th 
qubit: = —ijpaz, with 1^ being the absolute value of 
the persistent current in the z-th qubit; then the eigenstates of 
a z correspond to the clockwise (a z |t) = — |4-}) and counter- 
clockwise (er 2 |t) = |t)) current in the i-th qubit. The energy 
bias Ci(t) is controlled by constant and alternating magnetic 
fluxes 

ti(t) = 2/W U t (t) - i$ ) = 4 0) + Ut), (2a) 

€i(t) = 24 l )$ /acSin^, /ac = |p. (2c) 

The basis state vectors for the two-qubit system 
{|44)> I4t)> \ti}> Iff)} are composed from the single-qubit 
states: |tt) = |4-)(i)|t)(2)> etc. For identification of the level 
structure and understanding different transition rates, we will 
start the consideration from the case of two non-interacting 
qubits. Then, the energy levels of two qubits consist of the 
pair-wise summation of single-qubit levels, 

^ = ±^ = 4vM 0)2 + A?, (3) 
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FIG. 1: (Color online). Schematic diagram of the two-qubit sys- 
tem. Two different flux qubits are biased by independent constant 
magnetic fluxes, and , and by the same alternating mag- 
netic flux $ ac sin tot. The former controls the energy levels structure, 
while the latter changes the populations of the levels. The dissipa- 
tion processes are described by coupling the system to the bath of 
harmonic oscillators. 



which are the eigenstates of the single-qubit time-independent 
Hamiltonian (fl]i at / ac = 0. We demonstrate this in Fig. [2 a), 
where we plot the energy levels, fixing the bias in the first 
qubit fx, as a function of the partial bias in the second qubit 
/2- Then the single-qubit energy levels appear as (dashed) 

horizontal lines at Ef = ±\\J e^ 2 + Aj for the first qubit 

and as the parabolas at E% (/ 2 ) = ±\\J e 2 0) {h) 2 + A§. 

After showing the two-qubit energy levels in Fig. 12a), we 
assume that the relaxation in the first qubit is much faster than 
in the second (this will be studied in the next Section), which 
is shown with the arrows in the figure. And now our prob- 
lem, with four levels and with fast relaxation between certain 
levels, becomes similar to the one with lasers. 7 This allows us 
to propose three- and four-level lasing schemes in Fig.|2jb,c). 
This is the subject of our further detailed study. 

We have analyzed the relaxation in the system of two un- 
coupled qubits. However this system can not be used for las- 
ing, since this requires pumping from the ground state to the 
upper excited state (see Fig. |2jb,c)). Such excitation of the 
two-qubit system requires simultaneously changing the state 
of both qubits and can be done provided the two qubits are 
interacting. That is why in what follows we consider in de- 
tail the system of two coupled qubits. The coupling between 
the two qubits we assume to be determined by an Ising-type 
(inductive interaction) term ^a^a^ 2 \ where J is the cou- 
pling energy between the qubits. Then the Hamiltonian of the 
two driven flux qubits can be represented as the sum of time- 
independent and perturbation Hamiltonians 

H 2q = Ho + V(t), (4) 

H a =j2 - s^ ) + r^ 2) ' (5) 

i=l,2 ^ ' 
i=l,2 



2 




FIG. 2: (Color online). Energy level structure of two uncoupled 
qubits (J = 0). (a) One-qubit and two-qubits energy levels are 
shown by dashed and solid lines as a function of partial flux fa at 
fixed flux /i. We mark the energy levels by the current operator 
eigenstates, |44) etc. Particularly, we will consider the energy levels 
and dynamical behaviour of the system for the flux biases fa = fai^ 
(marked by the square) and fa = /2R (marked by the circle). By the 
arrows we show the fastest relaxation - for qubit 1. (b) Scheme for 
three-level lasing at fa = fah- The driving magnetic flux pumps (P) 
the upper level |3). Fast relaxation (R) creates the population inver- 
sion of the first excited level 1 1) in respect to the ground state |0); 
these two operating levels can be used for lasing (L). (c) Scheme for 
four-level lasing at fa = fan.. Pumping (P) and fast relaxations (Ri 
and R2) create the population inversion of the level |2) with respect 
to level |1). 

where ct^I = &x,z ®gq, a^l — ao ® a x , z , and <tq is the 
unit matrix. When presenting concrete results we will use the 
parameters of Ref.[§ Ai/ft = 15.8 GHz, A 2 /h = 3.5 GHz, 
I^^o/h = 375 GHz, I^^o/h = 700 GHz, J/h = 3.8 
GHz. 

For further analysis of the system, we have to convert to 
the basis of eigenstates of the unperturbed Hamiltonian (O. 
Eigenstates {|0), |1), |2), |3)} of the unperturbed Hamilto- 
nian ((5)1 are connected with the initial basis 
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where S is the unitary matrix consisting of eigenvectors of 
the unperturbed Hamiltonian Making use of the trans- 
formation Hq = S^^^HqS, we obtain the Hamiltonian H' in 
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the energy representation: H' =diag(i?o, E%, E2, £3). These 
eigenvalues of the Hamiltonian H are computed numerically 
and plotted in Fig. [3a) as functions of the bias flux in the 
second qubit f%. The distinction from Fig. 0a), calculated 
with J = 0, is in that, first, the crossing at = /| be- 
comes an avoided crossing, and second, the distance between 
the [previously single-qubit] energy levels is not equal, e.g. 
now E3 — E 2 ^ Ei — Eq. 
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FIG. 3: (Color online), (a) Energy levels of the system of two cou- 
pled qubits. Arrows show the pumping and dominant relaxation, as 
in Fig. |2] (b) The relaxation rates Wmn, which give the probability 
of the transition from level n to level m, induced by the interaction 
with the dissipative bath. Dominant relaxations are W13 and W02 to 
the left from the avoided crossing at f% = /J and W23 and Wqi to 
the right. (The small relaxation rates W03 and W12 are not shown.) 

Likewise, we could also convert the excitation operator 
V(t) to the energy representation 



v'(t) = s- i v(t)s= -hm 



(8) 



(9) 



III. MASTER EQUATION AND RELAXATION 

A. Bloch-Redfleld formalism 

Following Ref. 10, we will describe the dissipation in the 
open system of two qubits, assuming that it is interacting with 
the thermostat (bath), see Fig. Q] Within the Bloch-Redfield 
formalism, the Liouville equation for the quantum system in- 
teracting with the bath is transformed into the master equation 
for the reduced system's density matrix. This transformation 
is made with several reasonable assumptions: the interaction 
with the bath is weak (Born approximation); the bath is so 
large that the effect of the system on its state is ignored; the 
dynamics of the system depends on its state only at present 
(Markov approximation). Then the master equation for the 



reduced density matrix pit) of our driven system in the en- 
ergy representation can be written in the form of the following 
differential equations 10 



Pij = -mjPij - ^ 



V , p + 5jj } j p nn Wj n "/ij Pij ■ 

(10) 



Here tUij = (Ei — Ej ) /h 7 and the relaxation rates 

*^mn — 2 Re r n rn m n , 
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are defined by the relaxation tensor Ti mn k, which is given by 
the Golden Rule 

00 

Ti mnk = i J dte- l ^ kt (H Um (t)H link (p)) . (13) 



Here H\(t) is the Hamiltonian of the interaction of our sys- 
tem with the bath in the interaction representation; the angular 
brackets denote the thermal averaging of the bath degrees of 
freedom. 

It was shown 9 ' 11 that the noise from the electromagnetic 
circuitry can be described in terms of the impedance Z(u) 
from a bath of LC oscillators. For simplicity one assumes 
that both qubits are coupled to a common bath of oscillators, 
then the Hamiltonian of interaction is written as 



H, = 



(14) 



in terms of the collective bath coordinate X = J2k c k^k- 
Here stands for the magnetic flux (generalized coordinate) 
in the k-th oscillator, which is coupled with the strength Ck to 
the qubits. We note that the coupling to the environment in 
the form of Eq. ( TT4b applies only to correlated noise, or both 
qubits interacting with the same environment. One could ar- 
gue that it would be more realistic to use two separate terms, 
one for each qubit coupled to its own environment. However, 
since this term leads to different relaxation rates in our qubits 
1 and 2 (see below), then the form in Eq. (fT4l) should give 
essentially the same results as two separate coupling terms. 

Then it follows that the relaxation tensor Ti mn k is defined 
by the noise correlation function S{uS) 



S(u) = / dte- wt (X(t)X(0)) 



nk 



(15) 



(16) 



(17) 
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The noise correlator S (w) was calculated in Ref. 1 1 within 
the spin-boson model and it was shown that its imaginary part 
results only in a small renormalization of the energy levels 
and can be neglected. The relevant real part of the relaxation 



tensor 



Rel^,^ — — Ai mn k J (<jJ n k) 



coth ■ 



IT 



1 



(18) 



is defined by the environmental spectral density J(w). Here T 
is the bath temperature (ks is assumed 1); for the numerical 
calculations we take T/h = 1 GHz (T = 50 mK). The elec- 
tromagnetic environment can be described as an Ohmic resis- 
tive shunt across the junctions of the qubits, Z{uj) — R. 9 Then 
the low frequency spectral density is linear J(ui) oc ujZ(ui) oc 
uj and should be cut off at some large value w c ; the realistic 
experimental situation is described by 11 



J(«) 



1 + oj 2 /u) 2 



(19) 



where a is a dimensionless parameter that describes the 
strength of the dissipative effects; in numerical calculations 
we take a = 0.01 and cj c /27t = 10 4 GHz (the cut-off fre- 
quency lj c is taken much larger than other characteristic fre- 
quencies, so that for relevant values lu : J(ui) w ahaj). 



B. Relaxation rates 

From the above equations the expression for the relaxation 
rates from level \n) to level \m) follows 



W, 



- — A 
4ft 



huj mn 

coth 1 

2T 



(20) 



These relaxation rates are plotted in Fig. Ob) as functions 
of the partial flux bias /2. This figure demonstrates that the 
fastest transitions are those between the energy levels corre- 
sponding to changing the state of the first qubit and leaving 
the same state of the second qubit, cf. Fig. Oa). Namely, 
the fastest transitions are those with the rates W13 and W02 
to the left from the avoided crossing and W23 and Wo% to the 
right, which correspond to the transitions |tt) — > |4-t) and 
|t4-) — ► I44-)- Note that we do not show in the figure the rates 
W03 and Wx%; they correspond to the transitions with simul- 
taneously changing the states of the two qubits and they are 
much smaller than the rates shown. 

The relaxation rates Wij are shown in Fig. |4]as functions 
of the two partial bias fluxes, fx and fa. Again, one can see 
the regions where certain relaxation rates are dominant. Such 
a difference in the relaxation rates creates a sort of artificial 
selection rules for the transitions similar to the selection rules 
studied in Refs. 



1243. 



In our case the transitions are induced 
by the interaction with the environment and the difference is 
due to the different parameters of the two qubitsJi To further 
understand this issue, we consider the single-qubit relaxation 
rates. 

From the above equations we can obtain the energy relax- 
ation time T\ and the decoherence time T2 for single qubit. 
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FIG. 4: (Color online). Relaxation rates Wmn versus partial biases 
of the two qubits, f± and fa. The square and the circle show the 

at which the calculations of other 



parameters /1 and f 2 = / 2 l(r) 
figures are done. 



For the two-level system with two states |0) and |1) the relax- 
ation time is given by 10 — Wqi + W\q. The Boltzmann 
distribution, W10/W01 = exp(— AE/T), means that at low 
temperature the major effect of the bath is the relaxation from 
the upper level to the lower one. Now, from Eq. ( f20b it follows 
that 

, aA 2 AE 

7T 1 = — t— coth . (21) 

1 2HAE 2T 

Also from Eq. ( fT2l we obtain the dephasing rate 10 

1 olT e(°) 2 

T 2 = Re ^ = 2 T " + -Ai^ (22) 
For the calculation presented in Fig. |2|a) for two qubits with 
J = in the vicinity of the point = f% , where AE^ — 
AE( 2 \ we obtain 



(i) 



,(2) 



A2 

Ai 



(23) 



As we explained above, the lasing in the four-level system 
requires the hierarchy of the relaxation times. In particular, we 



assumed t[ 1] < T^' . So, in our calculations we have taken 
Ai ^> A2 and consequently the first qubit relaxed faster. This 
qualitatively explains the dominant relaxations in Fig.|3b). 



C. Equations for numerical calculations 

If we use the Hermiticity and normalization of the density 
matrix, then the 16 complex equations ( TTOb can be reduced to 
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15 real equations. After the straightforward parametrization 



of the density matrix, pij 

[V',y] 



l 

h 



w ir x r 

Vij = -UijXij + \ [V',y]i 



+ iy tj , we get 1J 



- my a, i> i; 



1,2,3; 

(24a) 
(24b) 

(24c) 



= 0, Xqo = 1 - (xn + x 2 2 + x 33 ); x jt = Xij, y^ = -y ir 
This system of equations can be simplified if the relaxation 
rates are taken at zero temperature, T = 0, and neglecting 
the impact of the inter-qubit interaction on relaxation, J = 0. 
Then among all the Wij and jij non-trivial are only the ele- 
ments corresponding to single-qubit relaxations (see Eqs. (EH 
l22l ). For example consider / 2 < / 2 (see Fig. |2ja) for the 
notation of the levels), then non-trivial elements are 



W 13 = W n2 = 



23 



W 



01 



(If))" 1 
(if))" 1 



2ftA£ 2 ' 



(25a) 
(25b) 



713 
723 



731 = 702 = 720 



Pr)- 1 = -(T 1 (1 V 1 (26a) 
732 = 7oi = 710 = (if)" 1 = i(T 1 (2) )- 1 (26b) 



In our numerical calculations we did not ignore the influ- 
ence of the coupling on relaxation, i.e. we did not assume 
J = 0. However, we have numerically checked that such sim- 
plification, J = 0, resulting in the relaxation rates (25-26), 
sometimes allows to describe qualitatively dynamics of the 
system. 



IV. SEVERAL SCHEMES FOR LASING 

In Sec. II and in Fig. [2] we pointed out that in the system 
of two coupled qubits there are two ways to realize lasing, 
making use of the three or four levels to create the popula- 
tion inversion between the operating levels. In this Section 
we will demonstrate the lasing in the two-qubit system solv- 
ing numerically the Bloch-type equations (l24b with the relax- 
ation rates given by Eqs. ( fTTl[T2l[T8] l. Besides demonstrating 
the population inversion between the operating levels, we ap- 
ply an additional signal with the frequency matching the dis- 
tance between the operating levels, to stimulate the transition 
from the upper operating level to the lower one. So, we will 
first consider the system driven by one monochromatic signal 
f(t) = f ac svnojt to pump the system to the upper level and 
to demonstrate the population inversion. Then we will apply 
another signal stimulating transitions between the operating 
laser levels: 



Solving the system of equations d24l . we obtain the population 
of z-th level of our two-qubit system, Pi = x«. The results of 
the calculations are plotted in Figs. [5] and [6] where the tempo- 
ral dynamics of the level populations is presented for different 
situations. 
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f(t) = /ac sinwt + / L sin Lo h t. 



(27) 



FIG. 5: (Color online). Three-level lasing and stimulated transi- 
tion. Time evolution of the numerically calculated occupation prob- 
abilities at biases fi = 14 X 1CT 3 and fa = 11 X 10~ 3 is plotted 
for (a) one-photon driving and (b) two-photon driving. As shown in 
the inset schemes, the driving and fast relaxation create the inverse 
population between the levels |1) and |0). So, these levels can be 
used for lasing, which we schematically mark by the double arrow. 
After some time delay (when the population inversion is reached) 
an additional periodic signal (S) /l cos LJht is turned on matching 
the operating levels, huJh — E\ — Eq. This leads to the stimulated 
transition |1) — > 0). 



In Fig. [5] we consider the situation where the relevant dy- 
namics includes three levels (for definiteness, we take fi — 
14 x 1CP 3 , /2 = 11 x 10~ 3 , which is marked as the square 
inFig.lUi. Pumping (| 0} ->■ |3)) and relaxation (|3) -J- |1}) 
create the population inversion between the levels |1) and |0). 
For pumping we consider two possibilities: one-photon driv- 
ing, Fig. |3 a), when hcu = E 3 — Eq, and two-photon driving, 
Fig-Hb), when 2hu> = E 3 — Eq. In the latter case we have 
chosen the parameters (namely /i and / 2 ) so, that the two- 
photon excitation goes via an intermediate level |2). We note 
here that, as was demonstrated in Ref. 8, the multi-photon 
excitation in our multi-level system can be direct, as below 
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FIG. 6: (Color online). Four-level lasing and stimulated tran- 
sition. Time evolution of the occupation probabilities at biases 
fi = 14 x 10" 3 and f 2 = 20 x 1CT 3 is plotted for (a) one-photon 
driving and (b) two-photon driving. The driving and fast relaxation 
create the inverse population between the levels |2) and |1). After 
a time delay an additional periodic signal /lCoswl£ is turned on 
matching the operating levels, /kjl = E% — E\. This leads to the 
stimulated transition |2) — ¥ |1). 



in Fig. |6jb), or ladder-type, via an intermediate level, as in 
Fig. |3b). Figure [5] was calculated for the following param- 
eters: wl/27t = 13.7 GHz (Swl = E\ — Eq) and also (a) 
a;/2vr = 35.2 GHz, / ac = 7 x 1CT 3 , f L = 5 x 1CT 3 ; (b) 
w/2vr = 17.6 GHz, / ac = 2 x 1CT 3 , f L = 5 x 1(T 3 . 

Next, we consider the scheme for the four-level lasing, 
which occurs in a similar scenario, except the changing of the 
levels. Then, the main relaxation transitions are |3) — > |2) 
and |1) — > |0), and now the population inversion should be 
created between levels |2) and |1). For this we take the par- 
tial biases fx = 14 x 1CT 3 , f 2 = 20 x 10~ 3 (marked by 
the circle in Fig. |4j. First, the system is pumped only with 
one signal either with Twj = E$ — Eq, Fig. |6ja), or with 
2hw = E3 — Eq, Fig. |6jb). Such pumping together with fast 
relaxation (|3) — > |2)) creates the population inversion be- 
tween the levels |2) and |1). Fast relaxation from lower laser 
level |1) into the ground state |0) helps creating the population 
inversion between the laser levels |2) and |1), which is the ad- 
vantage of the four-level scheme. 7 Then the second signal is 
applied with a frequency matching the laser operating levels 



(hoj-L = Ei — Ei). This stimulates the transition |2) — > |1), 
which provides the scheme for the four-level lasing. Figure 
[6] was calculated for the following parameters: = 9 

GHz (hu L = E 2 - E x ) and also (a) w/2tt = 47.4 GHz, 
/a C = 5 x 1(T 3 , / L = 3 x 1CT 3 ; (b) w/2vr = 23.7 GHz, 
/ ac - 5 x ltr 3 , / L = 5 x 10- 3 . 

In the experimental realization of the lasing schemes pro- 
posed here, the system of two qubits should be put in a quan- 
tum resonator, e.g. by coupling to a transmission line res- 
onator, as in Ref. 5. Then the stimulated transition between 
the operating states, which we have demonstrated here, will 
result in transmitting the energy from the qubits to the res- 
onator as photons. For this, the energy difference between 
the operating levels should be adjusted to the resonator's fre- 
quency. 



V. CONCLUSIONS AND DISCUSSION 

We have considered the dissipative dynamics of a system of 
two qubits. Assuming different qubits makes some of the re- 
laxation rates dominant. With these fast relaxation rates, pop- 
ulation inversion can be created involving three or four levels. 
The four-level situation is more advantageous for lasing since 
the population inversion between the operating levels can be 
created more easily. We demonstrated that the upper level can 
be pumped by one- or multi-photon excitations. We also have 
shown that after applying additional driving, the transition be- 
tween the operating levels is stimulated. 

When presenting concrete results, we have considered the 
system of two flux superconducting qubits with the realistic 
parameters of Ref. |8j For lasing in a generic two-qubit (four- 
level) system, our recipe is the following. The hierarchy of 
the relaxation times in the system is obtained by making it 
asymmetric, with different parameters for individual qubits. 
This makes transitions between the levels corresponding to a 
qubit with smaller tunneling amplitude A negligible, which 
creates a sort of the artificial selection rule. Based on our 
numerical analysis, we conclude that the optimal combination 
of pumping and relaxation is realized for Ai ^> A2 ~ J. 

Creation of the population inversion and the stimulated 
transitions between the laser operating levels, demonstrated 
here theoretically, can be the basis for the respective experi- 
ments similar to Ref. 5. In that work, a three-level qubit (ar- 
tificial atom) was coupled to a quantum (transmission line) 
resonator. First, spontaneous emission from the upper oper- 
ating level was demonstrated. In this way the qubit system 
can be used as a microwave photon source^ Then, the op- 
erating levels were driven with an additional frequency and 
the microwave amplification due to the stimulated emission 
was demonstrated. We believe that similar experiments can 
be done with the two-qubit system (which forms an artificial 
four-level molecule from two atoms/qubits). To summarize, 
we propose to put the two-qubit system in a quantum res- 
onator with the frequency adjusted with the operating levels 
and to measure the spontaneous and stimulated emission as 
the increase of the transmission coefficient. Such lasing in a 
two-qubit system may become a new useful tool in the qubit 
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